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Abstract One of the ASACUSA (Atomic Spectroscopy And Collisions Using 
Slow Antiprotons) collaboration’s goals is the measurement of the ground state 
hyperfine transition frequency in antihydrogen, the antimatter counterpart of 
one of the best known systems in physics. This high precision experiment yields 
a sensitive test of the fundamental symmetry of CPT. Numerical simulations 
of hyperfine transitions of antihydrogen atoms have been performed providing 
information on the required antihydrogen events and the achievable precision. 
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1 Introduction 

Through a precise measurement of the hyperfine splitting of antihydrogen and 
a comparison with hydrogen, the ASACUSA collaboration aims at testing the 
CPT symmetry pQ. In ground state hydrogen the interaction between proton 
and electron spin leads to a singlet state with quantum number F = 0 and a 
triplet state with F = 1 (see figure [I]). The transition frequency between these 
two levels is one of the most accurately measured quantities and is therefore 
well suited to test CPT with very high precision m- In the presence of a 
magnetic field the degenerated F = 1 level splits up (see figure [l]) and the 
energies of all four states shift. This is described by the Breit-Rabi formulae 
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Fig. 1 Is groundstate levels of antihydrogen and level splitting in a magnetic field, as well 
as the two transitions tt\ (F = 1 , Mp = — 1 ) —>• (F = 0, Mp = 0) and <j\ (F = 1 , Mp = 
0) ->► (F = 0 ,M f = 0). 


0 : 

£(i,i) = 4^0 ^ 2 ( gj + 9 i )^ bB (1) 

£ 7 (i,o) = ^4^0 + -^oVl + a-’ 2 (2) 

£(i,-i) = 4^0 ^ ^(gj + gi)g B B (3) 

£( 0 , 0 ) = -^£o - ^£o\/l +X 2 (4) 

where Eq = hv o with the zero magnetic field transition frequency bq and the 
shifted energies £( 1 , 1 ) etc. of the four hyperfine levels. B denotes the external 
magnetic field, x = B/Bq and B o = 2irvo/((gj — gi)g B ) with the electron 
g-factor gj = —2.0023193043718 [5] and <?/ = g p m e /m p = 0.003042064412 [5] 
with the proton g-factor g p (both in units of the Bohr magneton /is). Conse¬ 
quently not only one but several transitions can be observed. These states can 
be classified into low- and high-field seekers (LFS and HFS) depending on their 
behaviour in an inhomogeneous magnetic field. Depending on the alignment of 
their mangetic moment in a field, atoms with parallel magnentic moment will 
tend to move toward higher field regions, for states with antiparalell magnetic 
moment the opposite is the case. Relevant for the ASACUSA experiment are 
the fji and the transition: 

fTi : (1, 0) -)■ (0, 0) : u ai = v^\J\ + x 2 (5) 

TTi : (1,-1) -t (0,0) : v Wl = - \{gj + gi)\i B B/h + ^ 0 a/1 + x 2 (6) 

Numerical simulations of these transitions have been done in order to de¬ 
termine the obtainable precision and will be discussed in the next section. The 
experimental resolution of the ground state hyperfine transitions is inversely 
proportional to the interaction time of the atoms with the microwave field in 
the cavity. The line width of a resonance scan can therefore be estimated by 
5u = 0.799 /T |6j where T is the flight time of the atoms through the cavity. 
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The length of the cavity is 10 cm and the average beam velocity is expected to 
be 1000 m/s which leads to a estimated FWHM of 7x 10 -6 . With good enough 
statistics it should be possible to determine the center of the resonance spec¬ 
trum, which corresponds to the transition frequency, with a relative precision 
of 1CT 7 . 


2 Simulations 

The ASACUSA spectrometer line is built up of an antihydrogen source [7), 
a spin flip inducing microwave cavity, a spin analyzing sextupole magnet and 
a detector. The oscillating magnetic field B osc which induces the spin flips 
is provided by the double stripline resonator of the cavity placed inside a 
cylindrical vacuum tank [5j. B osc is very homogeneous in the plane orthogonal 
to the beam and has a sinusoidal distribution parallel to the beam [5j. As a 
consequence of this, the resonance spectra of a transition will have a double 
peak structure leading to zero signal (i.e. no induced spin flips) when the 
frequency of the oscillating field is on resonance. This can be seen in the plot 
of the simulated resonance scan as shown in figure [2] 

In the ASACUSA cavity two kinds of transitions are possible - the cri and 
the 7ri - depending on the angle between the oscillating magnetic field and the 
static magnetic field, B sta t , provided by Helmholtz coils. For the a± transition, 
the two B-fields need to be parallel whereas for the transition they have to 
be orthogonal to each other. In the experiment a small static magnetic field 
is used. Between 0 and 10 Gauss the transition frequency of oy has a second 
order dependence on the external field whereas the 7Ti transition varies linearly 
with the external field and is thus more sensitive to B-field inhomogeneities. 

There are now two possible ways of determining the transition frequency at 
zero static magnetic field B sta t • Firstly, by measuring the resonance frequency 
of the (7 1 transition at different external B-fields B\ tat and then extrapolating 
to zero field using the Breit-Rabi formula or secondly, by measuring a reso¬ 
nance scan of both the ay and 7iy transition at the same B stat and using the 
Breit-Rabi formulae in order to extract the hyperfine transition frequency: 



i/ 0 = 


where g± = gi ±gj. Simulations of the setup are being done using the particle 
physics toolkit Geant4. It is designed for high energy physics and therefore 
had to be modified in order to allow simulations of hyperfine transitions in 
a microwave field. The evolution of the spin state of the atoms for a certain 
time in the radiation field of the cavity is determinated by solving the op¬ 
tical Bloch equations using a Runge Kutta algorithm with adaptive stepsize 
control. One of the simulated resonance scans can be seen in figure [2] For rea¬ 
sons of simplicity two gaussians were used as a fit function to determine the 
center position which corresponds to the transition frequency. A comparison 
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Fig. 2 Resonance scan of the ry\ transition: count rate on the detector in dependance of 
microwave frequency. A beam with 3.5x 10® detected particles per scan, Maxwell-Boltzmann 
distributed velocity with a temperature of 50 K, a polarization of 70% HFS states and 30% 
LFS states and a measured fieldmap with B\ v = 8.78x10“ 5 T with rmsi = 5.362x10“' T. 


of this simple fit function with a more sophisticated fit algorithm showed no 
significant difference in precision and value for the center frequency, hence the 
simpler fit is sufficient for the goal of the present study. 

The two possible ways of determining the zero field transition frequency 
have been compared. When only using the oq transition, scans have been sim¬ 
ulated at three different static magnetic fields B l stat using measured fieldmaps 
inside the cavity, in order to make the simulations as realistic as possible. 
The mean field strength and RMS of the measured fieldmaps used are: B\ v = 
8.78x10“® T with rmsi = 5.362xl0“ 7 T, B 2 av = 2.194xlCT 4 T with rms 2 = 
1.414xl0“ 6 T and B% v = 4.389xl0“ 4 T with rms 3 = 2.911x10“® T. A beam 
with a Maxwell-Boltzmann distributed velocity of 50 K and a polarization of 
70% LFS and 30% HFS has been used, since this is the expected polarization 
of ASACUSA’s antihydrogen source (TO] . 

For the second method, the same beam properties were used. The angle be¬ 
tween B stat and B osc was changed to 45° in order to have a parallel and a 
orthogonal component to B osc and being able to drive both resonances in the 
same configuration. 

With very high statistics, (3.5x10® atoms per scan detected) as in figure [ 2 ] 
the method using only the ay transition showed a larger error by a factor of 
«1.12. When using lower statistics the method using only the oq transition 
becomes more favorable due to the sensitivity on inhomogeneities of the 
transition. Simulations with 2.1x10® particles per scan reaching the detector, 
which is the particle number where the sensitive 7Ti transition becomes visible, 
have been done. Comparing the relative error of the two methods shows that 
the second one yields a higher relative error by a factor of «1.4 for the homo- 
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Fig. 3 Extrapolation to 0 field using simulations at 3 different fields and 35 points per 
scan: Absolute (points) and relative errors (triangles) depending on particle number per 
scan. Two gaussians were used as a fit function. 


geneity of the previously listed measured field maps. This can be improved by 
decreasing the inhomogeniety of B stat which is already planned. 

In figure [3] the achieved precision of the zero magnetic field HFS is plotted as 
a function of the number of detected particles. Here, the method of extrapo¬ 
lating three <j\ transitions at different B\ tat was used. 

Simulations show that the lowest possible particle number at which the dou¬ 
ble peak structure is still visible is about 2800 detected atoms per scan when 
using 35 points yields a relative precision of « 4.5xl0 -7 . The relative error 
corresponding to 3.5xl0 6 detected atoms per scan improves this by an order 
of magnitude. 


3 Summary 

Results of numerical simulations of hyperfine transitions of antihydrogen within 
the ASACUSA cavity have been presented. The zero static magnetic field tran¬ 
sition frequency has been determined by two different methods using the a\ 
and the 7Ti transition frequencies and the resulting precisions have been com¬ 
pared. Simulations show that using both transitions only leads to a higher 
precision when using high statistics due to the sensitivity to magnetic field 
inhomogeneities of the n 1 transition. The impact of particle numbers on the 
determined hyperfine transition frequency and its relative error has been dis¬ 
cussed. Simulations show that the minimum particle number needed, in order 
to see the characteristic double peak structure, is 2800 atoms per scan and the 
achievable relative error ss 4.5 x!0~ 7 . 
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